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Recently the Casmir-Polder force felt by an atom near a substrate under nonequilibrium stationary 
conditions has been studied theoretically with macroscopic quantum electrodyanamics (MQED) and 
verified experimentally with cold atoms. We give a quantum field theory derivation of the Langevin 
equation describing the atom's motion based on the influence functional method valid for fully 
nonequilibrium (nonstationary) conditions. The noise associated with the quantum field derived 
from first principles is generally colored and nonlocal, which is at variance with the 'local source 
hypothesis' of MQED's generalization to nonequilibrium conditions. Precision measurements on the 
shape deformation of an atomic gas as a function of its distance from a mirror would provide a 
direct check of our predictions based on this Langevin equation. 
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The Casimir-Polder (CP) force [T] between a neutral 
atom and a mirror or a dielectric surface [2] has drawn 
renewed attention of theorists [3JH1IS] because of real pos- 
sibilities of detection [B|. At short distances (<100 nm) 
CP-Lifshitz type forces dominate the interaction between 
neutral bodies making them a relevant or even essential 
factor in the design of micromechanical devices, traps for 
cold atoms and in precision measurements for the detec- 
tion of deviations from known forces. 

At a finite temperature this force has two components 
originating from the vacuum and thermal fluctuations 
of the electromagnetic field, the latter is known as the 
Lifshitz force. Two major theoretical approaches are 
used: quantum field theory (QFT) under external condi- 
tions [7], which starts from microphysics based on QFT, 
and macroscopic quantum electrodynamics (MQED) [2], 
which is in the form of a linear response theory with an 
added stochastic source. MQED has been skillfully ap- 
plied [3] |3J [S] to the situation of a substrate at a different 
temperature from the field under nonequilibrium station- 
ary conditions. However one key assumption of MQED, 
that the fluctuations are local, has not been justified. It 
remains a challenge to bridge these two approaches, to 
remove or justify such an assumption, and to general- 
ize to fully nonequilibrium conditions for non-stationary 
systems. This is the aim of our research program. 

In this letter we lay out the basic structure of such 
a theory based on quantum open systems conceptual 
framework and the influence functional (IF) formalism 
8 . Our nonequilibrium formulation recovers the well- 
known CP force on an atom when the trajectory of the 
atom is stationary [Eq. ( 10|11 )], and the thermal CP 
force in a finite temperature field [Eq. ( 12|13 )]. More 
importantly, it gives a first-principles derivation of a 
Langevin equation which describes the atom's stochas- 
tic motion. The noise in this equation is generally col- 
ored and nonlocal, at variance with the main assumption 
of MQED. Our result for the dispersion of an atomic 



cloud could motivate experiments designed to measure 
its change in shape as a function of its distance from the 
mirror. 

The atom and its trajectory in a quantum field with 
boundary We model the internal degrees of freedom (idf ) 
Q of an atom by a (3-dimensional) harmonic oscilla- 
tor with natural frequency f2. The atom moves on a 
trajectory z(t) in a quantum field A^, the electromag- 
netic vector potential, in the presence of a mirror. The 
dynamics of the system is determined self-consistently 
by allowing these three variables (Q, A^,z) to negotiate 
amongst themselves. Even for a stationary atom it is 
necessary in the set up of the problem to assume its po- 
sition z to be a dynamical variable so its resultant trajec- 
tory comes from the mutual interactions with the other 
two variables. The action describing the entire system 
is S[Q, z] = S Q [Q] + S E [Av] + S z [z] + S int [Q, A», A 
(subscript E stands for the electromagnetic field which 
serves as an environment) with the action for the oscilla- 
tor given by Sq[Q] = f / d\[Q(\) 2 - n 2 Q(X) 2 } where 
m is the oscillator's reduced mass and A parameter- 
izes its worldline. The photon field action is given by 
S E [A»] = -i J d A xF^ u F^ where F av = d )1 A u — d v A u is 
the field strength tensor. The action for the motion of the 

atom's center of mass M is Sz [z\ — J dX[^AIz (A) — V[z\] 
where V[z\ is an external potential. 

In the dipole approximation, the Hamiltonian for an 
atom interacting with the photon field takes the form 
—d-E where d is the dipole moment of the atom and E 
is the electric field. In this spirit we define the interac- 
tion action S int [Q, z, A n ] = qj dXQ k (X)E k [z ,J, (X)} where 
Greek indices denote spacetime components, for time, 
Roman indices will be reserved for purely spatial com- 
ponents, and the Einstein summation convention is used 
throughout. 

World Line Influence Functional Assume that at time 
ti n the quantum statistical state of the oscillator, tra- 
jectory and field is described by a density operator 
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p(U n ). This state is unitarily evolved from the initial 
time ti n to a later time tf > ij n , and can be expressed 
in terms of path integrals by considering matrix ele- 
ments in an appropriate basis. The overall influence of 
the field on the dynamics of the atom is obtained by 
coarse-graining over the field variables resulting in the 
field-reduced density matrix [5], p r (Q f,Q'p Zf, Zf',tf) — 

I dQ m dQ' m J dz in dz\ n VQ /<| VQ> f£f4 v *>* 

x pz(zin, ^i n ', tin)-7 7 [J M_ ) J v+ ], where Vk is the measure 
for a path integral over the space of functions. This in- 
troduces the influence functional (IF) T\J^ , J u+ ] [5]. 

For the coupling given above and assuming an initially 
uncorrelated and Gaussian state the influence functional 
can be calculated exactly and is given by 

x J dV[^(y,y')^ + (2/') + ^^(y,y')^(2/')]}-(i) 

Here the current density is J^{x) = — q J cZA(9o^> + 
d j m^)S 4 (x t " - z^(X))Q j (X), J+ = (J+ J')/2 and J- = 
J — J' are its difference and semi-sum, respectively, and 
Vfiv = diag{ — \, 1,1,1) is the metric for Minkowski space. 
D r u e J(y,y') and D^ v (y,y') are the retarded Green's func- 
tion and Hadamard function for the field, repectively. 
They can be expressed in the Feynman gauge in terms of 
the corresponding Green's function for a massless scalar 
field in Minkowski space as D r **{x,x') — rj^Gre^x, x') 
and D^(x,x') = <q^ v G H {x,x'). 

To find the combined influence that the oscillator 
and the field have on the trajectory we need to fur- 
ther coarse grain the oscillator degrees of freedom result- 
ing in the oscillator-reduced influence functional (ORIF), 
T z [z-. / z+]. 

r -> -> rQf&f _ _ 

T z \z~^\ = / dQ f dQ in dQ' m / VQVQ 1 

xe *(s Q [Q}-s Q lQ']) pQ{ Q in , Q' m ■ t m )T\J^ , r+] (2) 

Noting that we cannot trace over the oscillator 
variables in ^ explicitly for arbitrary field boundary 
conditions, such as in the presence of a mirror, we proceed 
via a perturbative expansion in powers of the coupling. 
Writing ([2| in a more suggestive form Tz[z + ,z^] = 
exp{iS in/ [z»+ , *m- ; - % , - % } /o [J+ , J-] I . ± =Q , 

which defines the influence action, 
S mf [z^+,z^;Qj,Q+] = -i\nT[J^+,r-], and the IF 

for a three dimensional harmonic oscillator, fo[j + ,j~]- 
To factor the exponent out of the path integral [Q k± (X)] n 

is replaced with fo\j + , 3~] |j±=o- For a 

Gaussian initial state j~] can be evaluated exactly 



M1+J-] = Mexp{iJ d\d\'[f-(\) ■ j+(\')g ret (X,\') + 
|i-(A) •i-(V)to(A, A')]} where 5ret (A, A') and 5ff (A, A') 
are the retarded and Hadamard Green's functions for 
a one dimensional harmonic oscillator with natural 
frequency f2, J\f is a normalization constant, and the 
dot product is taken with respect to a 3 dimensional 
Euclidean metric. 

Expanding Q to lowest order in the coupling 
and partially resumming gives the ORIF to sec- 
ond order in the coupling —ih\Tz[z^,z~] ~ 
S lnf [z»+ fg n-.-i^,-i£]f o g+J-]\ J±=0 allowing 

us to write the reduced density matrix describing the 
center of mass motion. 

p r {zf,zj;t f )= I dz+ n dz~ n I Vz+Vz 

xe^- s ^pz(z^ n ,zr n ,t m )Tz[z^,n (3) 

Atom's mean trajectory The complex norm of the 
ORIF, \p r \ oc exp{-/ dXdX'z k -(X)N kj (X, X')z j -(X')} is 
non-vanishing and strongly suppressed for large values 
of the off diagonal elements, z~ — z — z 1 , as is indica- 
tive of decoherence of the quantum trajectory. N k j is a 
symmetric positive definite kernel quantifying the noise 
in the oscillator and field. 

Decoherence of the system due to its interactions with 
the quantum fluctuations of the environment and oscil- 
lator permits the existence of a semi-classical limit for 
the oscillator's path through space. Using a saddle- 
point approximation to evaluate Q about its clas- 
sical solution, z^i(X) = z k , one can show that the 
semi-classical dynamics is determined from the varia- 
tion SScGEA[z k+ , z k ~\l&z3~(r)\ z k- =a = where the 
so-called coarse grained effective action is given by 
ScGEA[z k+ ,z k -] = S z [z] - S z [z*} -i\nFz[z+,r\. 

Varying Sqgea with respect to z~ we obtain the mean 
(semi-classical) equation of motion [TU] 

Mz k (r) + d k V[z(T)} = / fe (r) (4) 

where the effective force, f k (r) (including back-action ef- 
fects), has the form 

fk(r) = ^ dxj^ dX'S ij 8(X ~ T)d k K%4 

L H {X,X')D^{X),z a {X')) 

+g ret (X,X')D^(z a (X),z a (X'))^ (5) 

where Kj = dof]j + djTjQ. Take caution to evaluate the 
derivatives before the particle trajectory is placed into 
the various kernels. 

The influence or back-action force on the oscillator tra- 
jectory describes dissipation and radiation reaction as 
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well as the forces due to constraints on the field. The 
first two effects must be taken into account when atom 
motion comes into play. In the following we assume an 
appropriate form for V[z\ so that Q admits static solu- 
tions where dissipative effects may be ignored. 

Casimir — Polder Force The placement of a mirror in 
the 2 = plane constrains the transverse components of 
the electric field to vanish there, and will lead to forces 
on the atom. This boundary condition can be accom- 
modated by appealing to the method of images. Thus, 
a dipole near a mirror will be attracted to its image on 
the other side, a classical electrostatic treatment for a 
permanent dipole gives al/z 4 dependence. When finite 
light propagation time and quantum fluctuations are ac- 
counted for this attractive force takes a modified form, 
1/z 5 , in the far field limit where the distance from the 
mirror is much greater than the period of the oscillator 
(c=l). This is the Casimir- Polder force. 

From the Green's function point of view the field con- 
straint can be satisfied by pairing every Green's func- 
tion with an image term i.e. G(a) — > G(a) — G(a) 
where a(x, x') is Synge's worldfunction defined to be half 
the geodesic distance between x and x 1 and a{x, x') = 
a(x,x') + 2zz'. The new terms, Fj? p , due to the pres- 
ence of a mirror are responsible for the CP effect. 



F k CP 



(r) 



<f_ 
2 



— Ai 



ds 5 ij d k K?rf, 



g H (s)D^(a[z a (r),z a (r~ S )}) 



+g ret ( S )D^(a[z a (r),z a (T^ S )}) 



(6) 



where a prime denotes differentiation with respect to 
the second argument. To accommodate the boundary 
conditions on the field the tensor structure of the im- 
age term changes D a p{a) = —(r/ a p — 2z a zp)G{<7) where 
z a = (0,0,0,1). 

The effective force from the image takes on the general 
form 



\.l\J 



2mO 



ds 



da 



d 

da 



2-i 



cosfis G re t{a) +sinf2s Gh{&) 



(7) 



where d/da operates only on the Green's functions for 
the field and &k = d^a. To find an explicit expression 
for the Casimir-Polder force we evaluate ^ for a static 
trajectory, z M (r) = (r, z) and i M (r) = (1,0). We find an 
analytic expression for the CP force in the long time limit 
when the field has dressed the atomic ground state and 
the field and oscillator were initially in the their respec- 
tive ground states. The CP-force has two contributions 



F, 



p CP2 _ <i 

k Z2n 2 mSlz^ kz 



8f2z + 6(l -2ti 2 z 2 )f(2Slz) 



4fiz(2fiV - 3)g(2Qz) 



(9) 



where f(x)(g(x)) is the auxiliary function for the cosine 
(sine) integral function, and {A, B} is the anticommuta- 
tor of A and B. 

Here F CP1 is derived from the term containing the re- 
tarded Green's function for the field and so is responsible 
for the electrostatic contribution to the CP force and the 
near field behavior. F CP2 is the dispersive part of the 
force because it contains the field Hadamard function. 
For interactions linear in the oscillator coordinate as as- 
sumed here, in perturbation theory the quantum ampli- 
tude to go from the ground state to any but the first ex- 
cited state vanishes. Thus the agreement of the present 
HO results with previous results for two level atoms using 
energy gradient methods [TT] is not surprising. 

In the near and far field limits we recover the asymp- 
totic expressions 



Qz « i f: 



CP 



3q 2 



Clz » 1 



7 CP 



32Trmttz 4 

3g 2 
8TT 2 mn 2 z 5 



(10) 
(11) 



where our results agree with the literature if we identify 
the static polarizability, a, with g 2 /47rmf2 2 , this form for 
a can be argued by examining the static solutions to the 
classical equations of motion. 

Thermal CP force The form of the CP force in a ther- 
mal field can be taken directly from ^ with all Green's 
functions replaced with their appropriate finite temper- 
ature version. The assumption of an initially factorized 
density matrix allows us to independently choose the ini- 
tial oscillator and field state. Choosing the oscillator 
and field to be in thermal states of different tempera- 
ture (with inverse temperatures /?, j3 respectively) gives 
rise to two distinct thermal contributions to the CP force. 

The retarded Green's functions appearing in (|6) will 
not contribute to the thermal effects as they are state 
independent. Modifications due to an initially thermal 
state will arise from the Hadamard functions only. The 
thermal Hadamard function for the field can be found by 
imposing a periodicity condition on the imaginary time [? 
]. For a harmonic oscillator it can be calculated directly 
g ^( T , t') = coth(/3fi/2) cosfi(r - r')/m!l This gives 



F° m = coth(/3n/2)^ 



CPl 



(12) 



4 




FIG. 1: Plot of thermal C-P force in units of hq 2 /m, for 
0, = 1, in the long time limit against perpendicular distance z 
(in units c/Q) of the atom from a plane mirror and absolute 
temperature T (in (h£l/kB)K) 
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2m0 



x rT—Xi 

V P.V. / ds 
Jo 



k=—oo 



xsinfis dkd v d u iGi{{t + ik(3{z),z,z') 



(13) 



where z — (i, x, y, —z). We have included the generalized 
case of a field state of spatially nonuniform temperature 
i.e. f3 — ► j3(x) as it is is Gaussian in field variables. 

In the high temperature, long time and far field limit 
we arrive at 



F, 



CP/32 



3q 2 Vkz i 

167r/3mfi 2 z 4 



(14) 



Stochastic trajectory The IF can produce a Langcvin 
equation for the trajectory with deviations from the mean 
caused by the quantum field fluctuations. It is given by 



Ml k + d a d k V[ Zl ]S a - d a f k [ Zi ]z a = Zk[zi 



(15) 



The key link in identifying a classical stochastic source 
(noise) from a quantum field is provided by the Feynman- 
Vernon identity for Gaussian integrals [5] . The two-point 
function for this classical stochastic source is related to 
the noise kernel by 

({6[ Z Q (A)],0[^(A')]}) = p mn 9 H (X,\') 

xd k d jlK ^D^[z a (X),z a (X% (16) 

The Langevin equation enables us to calculate the dis- 
persion of the atom's trajectory {^Az^t)), which is de- 
fined as the effective distance from the mean value that 
an ensemble of stochastic realizations takes. As the noise 
kernel contains the Hadamard function for the field it is 
sensitive to the boundary condition at z = 0. The im- 
age term present will make the distribution of noise vary 
with the distance from the mirror and in turn the dis- 
persion in the atoms' positions as well. This manifests 



as a fractional change in volume of a gas of noninteract- 
ing atoms. If we trap the atoms in a harmonic poten- 
tial with frequency il k in the kth direction, such that 



Ifi 



fl | >> q jmQrMz then the dissipation can be 



ignored in the final expression for the dispersion and we 
can directly compute the far-field long-time limit. The 
dispersion in the z-direction is given by 



15g 2 



1 



1 



167r 2 mfW 2 (fj 2 - tt 2 ) 2 



(17) 



where f2 z is the trapping potential frequency in the pres- 
ence of a mirror. The parallel components can be ob- 
tained from (17) by dividing by —15 and substituting the 



trap potential frequency for the unperturbed dimension. 
The expression for the dispersion shows that the presence 
of the mirror leads to a focusing in the perpendicular di- 
rection and a broadening in the parallel directions. 

In conclusion we have derived from first principles the 
semiclassical and stochastic equations for an atom's mo- 
tion near a mirror under fully nonequilibrium conditions. 
As a quantum field theory derivation of the stochastic 
source is given there is no need and no room for a 'local 
source hypothesis' which generalizes MQED to nonequi- 
librium conditions. Precision measurements (see e.g., 
|12j ) in the shape deformation of an atomic gas near a 
mirror as a function of atom-mirror spacings would pro- 
vide a direct check against our theoretical predictions. 
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